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The paper presents results of theoretical studies of sorption processes in 
submicroporous and microporous elastic materials. The sorption is viewed 
as consisting of adsorption and absorption. Mathematical formulas 
describing the processes are derived in a unified way by using 
a thermodynamic approach. They handle, both absorption-like phenomena 
in submicropores and constrained multilayer adsorption in micropores. 
The latter is described by a new BET-type equation (LBET model). 
Applications of the resultant multiple sorption model (MSM) are 
discussed. The proposed models are useful tools for interpretation of 
empirical sorption isotherms in terms of structural and energetic 
parameters of sorbent-sorbate systems. 


1. INTRODUCTION 


A number of important sorbents (polymers in glassy state, crystalline 
polymers, natural coal), are porous materials exhibiting elastic properties [1]. 
A characteristic feature of their structure is a noticeable fraction of pores which 
are of molecular size (sub- and micropores) [2-6]. The sorption process in such 
materials consists in: penetration of sorbate molecules into elastic phase 
(absorption); adsorption of the molecules on a surface of large pores; location 
of sorbate in submicropores (intermediate absorptive-adsorptive processes). 

The mathematical description of equilibrium properties of such absorptive- 
-adsorptive systems was derived in ref.[7], in a thermodynamic way. The 
treatment follows the line of the theory of polymer solutions (Flory [8]), 


Mathematical modelling of adsorptive-absorptive processes 37 


however far generalizations of this theory was necessary, to take into account 
different size of particles of the sorbent components (polymer chains and quasi 
crystalline domains). 

A simplified approach to this problem is presented in refs.[7, 9], where only 
adsorption and absorption subprocesses are considered, with entropy modelling 
being stressed (a so-called dual-sorption model - DSM). A generalization of the 
polymer solution theory was carried-out, aimed at modelling the sorption 
equilibria in porous materials exhibiting elastic properties. Formal assumptions 
required for the theoretical analysis of the thermodynamics of sorbent-sorbate 
systems were collected, and general formulas for enthalpy and entropy changes 
in the sorption process (respecting peculiarities of the sorbent structure) are 
presented. This model was found to be adequate for sorption of small molecules 
like water and methanol [10]. However, its examination showed that absorption 
capacity is significantly lower than might be expected from chemical and 
structural properties of hard coals [10, 11, 12]. 

To give a formal explanation of such effects and obtain a model satisfactory 
for larger sorbate molecules, we elaborated an extended mathematical 
description of sorption equilibrium in submicroporous elastic materials (like 
glass and crystalline polymers, hard coal), including absorption and a spectrum 
of adsorption-like subprocesses. It is referred to as the multiple sorption model 
— MSM. Basic assumptions, formulas and properties of MSM (related to those 
of DSM) were discussed in refs.[13, 14]. A hypothetical elastic sorbent being 
across-linked chain polymer with randomly distributed sphere-like holes 
(submicropores) was considered. Effects of submicropore size distribution on 
the equilibrium sorption capacity were analyzed in details from the energetic 
side. A distribution function of the hole radius was assumed, and sorption 
subsystems were defined, each containing holes which require the same energy 
for the sorption. The model completed with the full entropy terms, respecting 
the sorbent structure heterogeneity (like in DSM), was presented 
in refs.[15, 16]. 

In order to make the approach applicable to gaseous sorptives we elaborated 
a new state equation enabling us to evaluate thermodynamic parameters of 
sorbates in sorption systems at near and over-critical temperatures [16, 17]. 
Further development of the model was directed to improvements of the 
modelling of the adsorption process with respect to constrained multilayer 
adsorption [16, 18, 19]. In refs.[16, 18] we showed that MSM can be relatively 
easy completed with a BET-like formula describing the adsorption process in 
porous materials with dominant fraction of small micropores (LBET model). 
The full version of MSM was used in the paper [19] to examine sorption 
properties of hard coals with methane and carbon dioxide. Heterogeneity of 
coal structure and multilayer adsorption were taken into account. Applicability 
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of the model to pure adsorption process was also shown (adsorption of CH, and 
CO, on activated carbons). The paper [20] develops the theory with certain 
non-BET effects being taken into account, thus enlarging the model 
applicability up to the highest pressure range. Possible ways for handling the 
microporous surface heterogeneity are outlined in refs.[16, 21]. 

The main idea of our approach to the sorption modelling is to find links 
between measurable sorption data and physical properties of the sorbent and 
sorbate viewed separately. These properties, represented by appropriate 
parameters of DSM and MSM, are assumed to be the same in any sorption 
system. It enables us to get more information on sorption mechanisms and 
sorbent structure, by examination of sorption properties of a particular sorbent 
with different probing sorbates. 

This paper presents the multiple sorption model in its recent form, including 
certain newest concepts. Possible applications of the theory and exemplary 
results are shown. 


2. MULTIPLE SORPTION MODEL — THEORETICAL BASIS 


Let us consider an equilibrium state of the sorption system corresponding to 
a pressure P and temperature T. It may be viewed as an effect of the two stage 
procedure (see Figure1): 

a) removing of sorptive molecules from their volatile state (P, T, V,) into 
a liquid-like reference state (P, T, V,), with appropriate pressure P, and 
molar volume V, equal to that occupied by the sorbate molecules in the 
sorption system 

b) isobaric, isochoric and isothermal mixing of the molecules with the pure 
adsorbent components 

Thus, the mathematical description of the sorption may be derived by 

considering free energy changes at the above stages. In particular, the polymer 

solution theory [8] may be applied to describe the mixing stage. 

In elastic-porous material the mixing process is significantly dependent on 
amount of micropores and submicropores [2, 3, 5, 7, 12-14]. In the theoretical 
treatment such pores may be viewed as irregular holes of molecular and 
submolecular size, randomly distributed within the elastic phase [13, 14]. 
Meso- and macropores affect sorption properties to a much smaller extent [2, 3] 
as their surface area is practically negligible while compared to that of 
submicro- and micropores. In such materials the sorption process may proceed 
by penetration of sorbate molecules among elastic chains (pure absorption), and 
by placing of molecules in holes (pores) properly enlarged to the molecule size. 
In the latter case we may distinguish larger holes (micropores), which may hold 
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one or more sorbate molecules with no changes in size. The corresponding 
process will be referred to as pure adsorption, while that in smaller holes 
(submicropores) is of an adsorptive-absorptive nature. Differences in original 
hole size (and in corresponding expansion ratios) imply energetic heterogeneity 
of the absorption-like processes. In turn, the adsorption heterogeneity is due to 
diversity of geometric and chemical properties of particular sites placed on the 
surface of micropores. 
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Fig.1. Virtual creation of the sorption system 


Let us consider the sorption system as consisting of a number of subsystems 
a, each containing m,, moles of sorbate molecules with the same molar energy 
Q, and the molar volume V, [14]. In the equilibrium state of such a system the 
following general formula is valid [7, 8]: 


RT «inl LEO -E m, = AH -TAS (1) 
JEST)’ 
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where AH and AS denote total enthalpy and entropy change due to the mixing; 
R — gas constant, P, T) — sorbate fugacity. In case of vaporous sorptives 
(T<<T,; T, — critical temperature of the sorptive) one takes f(P,T)=P, 
(Ps T)=Po, U=Uzg and V=Vm (Vm Po, Us — molar volume of sorbate in liquid 
state, its saturated vapor pressure and evaporation energy, respectively, at 
temperature 7). For sorption at near and above critical temperatures (e.g. 
methane and carbon dioxide in room temperature) the quantities P,(V,), f(P,T) 
and cohesion energy U,(V,) are proposed to be determined by using the fluid 
state equation derived in the paper [17], with V, treated as a fitting parameter of 
the model. 
The enthalpy change AH may by expressed in the following concise form 


AH => m,,-Q, + 6H, (2) 


where Q, denotes the main component of the molar sorption energy contributed 
to the system by a-th subprocess, and 6H, is a component of relatively small 
magnitude representing effects of sorbate-sorbate interactions and of sorbent 
surface expansion (see ref.[14]). 

In our approach the energy Q, is viewed as the effect of combination of 
cohesion energies of the sorbent and sorbate according to the Berthelot rule. It 
may be expressed as [7, 13, 14]: 


Q, =V, p: |w,8? -¢, -28,6,]+82} 8, = JU,/V, (3) 


Qe — fraction of surface of sorbent molecules in total surface of dry sorbent 
and pores; 

ô. 5,— solubility parameter of sorbent and of sorbate, respectively; 

U, — molar cohesion energy of pure sorbate with the molar volume V,; 

©, — surface expansion ratio of holes attributed to a-th subsystem 


2 
o.-1-[%] for R,,<R,; @,=0 forR,, 2R, (4) 


P 


Ra, R, — the radius of a-th type hole and of sorbate molecule, respectively 

C, — a factor correcting the adhesion energy due to non perfect sorbent- 
-sorbate contacts (geometric factor) and effects of specific interactions 
(chemical factor): 
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1 for R, =0 (pure absorption) 
TOONE for 0<R,<R, 
=R „a/R, J= - 
6 =0R,/R,) Z, +(2, C, -2,) eE forR, <R,, <R, 
p B 
Ly iC =% for R „ 2 (pure adsorption) 
(5) 


where: Zg is an empirical parameter determining a fraction of effective sorbent 
sorbate contact area (like that occuring among molecules absorbed in polymers) 
for the given hole radius Rg; 

Z, — value of C, averaged over holes larger than sorbate molecule; C, — polarity 
factor expressing effects of specific interactions of polar sorbates with active 
groups on the surface. 

For nonpolar sorbates C,=1, so as ¢, is equal to the fraction of the hole wall 
surface viewed as being in full contact with the sorbate molecule (like in 
polymer solutions), while the remaining part of the surface is treated as being in 
contact with vacuum. 

For the pure adsorption( pointed to with the subscript A) (Ra 2 R,) we have 
@,=0; and 
Q, = (6 =24 sc -28,5,9° ): vV, (6) 
As it can be seen, the adsorption is mostly probable on local niches and cavities 
where Z, is appropriately large, and — in case of polar sorbates — on active 
groups, where C, is high. According to (6) the surface magnitude affects the 
adsorption capacity in two opposite directions. The larger surface, the greater 
number of adsorption sites, but they become less attractive due to decrease 
ing. 

In turn, energetic interactions in pure absorption may be characterized by the 
Flory-Huggins parameter x,,, calculated by using the equation (1) with ¢,=1, 
@,=1 and 9=1 (as for polymers [8]): 


Lye = = (82 =28,8, +8?)= = derti (7) 


It should be emphasized that sorption energy in the above model is 
expressed by cohesion energies of sorbate (U, or 6,) and of sorbent (6,) — see 
the formula (3) for Q, — those are treated as being independent on the particular 
sorption system (ð, can be taken from the van Krevelen formulas [1]). 
Moreover, in case of gaseous sorbates, U, and V, are strictly intercorrelated by 
the pure sorbate state equation (ref.[17]). For the pure adsorption and 
absorptive-adsorptive subprocesses the parameter Q, is corrected by the fitting 
parameters Z,, Rg, Zz (with respect to constraints: Z,<1, Zg <1), but for the pure 
absorption ¢, =1 is assumed (see Eq.5), so the absorption capacity is closely 
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correlated with that of remaining subprocesses. This specific feature of the 
model tightens considerably an area of acceptable values for its parameters, 
which allows us to draw more reliable conclusions on sorption mechanisms in 
submicroporous materials [19]. 

In classical adsorption models Q, is often treated as an empirical parameter 
characterizing the sorption system. The Eq.(3) shows that this quantity is very 
sensitive to changes in the surface properties (represented by @.° and Z4), hence 
Q, is not appropriate to characterize the system, if prediction of sorption 
properties is required. The proposed model seems to be more appropriate for 
such purposes. 

The entropy term AS may be derived with a modified Flory method [7, 8, 10, 
15, 16], based on a lattice model of the mixture containing empty holes. 
Finally, the formula for the sorption isotherm in a-th subsystem attributed to the 
holes of Rj. size may be derived by differentiation (with respect to mpa) of the 
formula (1) including equations for AH and AS. 

In this paper the recent version of MSM is presented, respecting the 
heterogeneous structure of the sorbent. In particular, effects of quasi-crystalline 
domains on both AS and AH are taken into account as depending on their 
volume and surface fraction in the system. 

The isotherm equation has the form: 


vi tv, (l -+ 
int = w, fin, -h-i In(l—v, :Bu)-2 Ti odtattalcal], 


Bau Z l-vy +4 "aa 


pth: 1o| (prem j aig) 
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where indices a, b point to different subsystems, and 
IT = f/f, = P/P, denotes the relative fugacity of the sorbate (relative pressure for 
vapours), 
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Vpa» Vpb — Volume fraction of sorbate molecules sorbed in a-th, b-th subsystem, 
respectively 

Vha — Volume fraction of original size holes attributed to a-th subsystem; v4, Ve; , 
Ve2 — volume fraction of quasi-crystalline domains, crosslinked chains and non- 
crosslinked chains, respectively, in sorption system; Waw, — volume expansion 
ratio of holes attributed to a-th, b-th subsystem: 


3 
w, -1-[% for R,,<R,; w,=0 forR,,2R, 
P 

üe, — volume fraction of crosslinked chains in sorbent; xy, Xe: , Xe2 -average 
size of quasi-crystalline domains, crosslinked chains and non-crosslinked 
chains, respectively, related to volume of sorbate molecule; B,, — geometrical 
parameter representing an effect of shape of paired {d,d} quasi-crystalline 
domains on the configurational entropy, calculated as: 


Ba ee elas ©) 


where Z,~ characterizes a compactness of quasi-crystalline domains 
(Zwa€ [1,z/2]) — see Eq.(18) given in paper [10]; z — lattice coordination number; 
aq — “specific surface area” of quasi-crystalline domains related to that of 
sorbate molecule; 
2 1 ) 
a, Slem 2 AA (10) 
z 2 
L., — number of typical crosslinkages per crosslinked chain; Q, Qa — the main 
energetic parameter defined in Eq.(3) for subsystems a and b, Qp», Ọps — surface 
fraction of sorbate molecules placed in b-th, B-th, subsystem of holes, 
respectively; @, , @. — surface fraction of holes, elastic chains, respectively, in 
the sorption system; @,’ — surface fraction of elastic chains in dry sorbent; 
9, — surface fraction of quasi-crystalline domains in dry sorbent; 
@, @, — surface expansion ratio of holes of a-th, b-th subsystem (see Eq.1); 
Ew, Eba , Eg — energetic constant defined for each pair of molecules {a,b}, 
{b,a} , {B,b} by equation (15) in ref. [7]; 


Re: sR) 
Bam 04455, h-o {1 Becta), liu, an 
P 
E,, & are expressed by equation (2) with C,=1; A, — total surface area of 
pores related to that of the sorbate molecule (in moles); dy), — changes in 


free surface area of coal due to swelling of coal matter, computed like for 
isotropic expansion of spheres, i.e.: 
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BA, = 214, -Z0 -0,) mp E, g (12) 
b bY, 


and the derivative of A, in Eq.6 is: 
ddA, 2 w 2 Va W 
= [o - 20-0, )oué, |-20-0,%, 52 (13) 


am, 3 9, 

m, — total sorption capacity; Mpa — amount of sorbate in the a-th subsystem of 

a sorption system (in moles); v, — volume fraction of coal molecules in the 

sorption system. 

In order to evaluate the surface and volume fractions (@ and v) of the 
system components, the model has to be completed with a quantitative 
description of the pore structure. The submicropores can be treated as sphere- 
-like holes of left-truncated normal size distribution with the mean value Rav 
minimal radius Ramin and dispersion Gg, [13]. Thus, if the total volume of the 
submicropores V,,,, is known, their specific surface area a, may be evaluated. 
It was stated in our earlier research that a, is usually larger then that of water 
molecule (refs [6, 7, 9-12]). The micropores are of more irregular shape, 
therefore there is no simple and reliable relation between their volume and 
surface area. However, the following assumptions seem to be acceptable: 

a) placing of a sorbate molecule at a site on a micropore surface disables 
a part of the micropore space for further independent first layer adsorption, 
and the volume excluded in this way depends on the sorbate molecule 
volume (micropore volume filling, rather than surface coverage); 

b) the surface of micropores may be roughly calculated by treating the first 
layer adsorption space (described above) as consisting of randomly placed 
spheres of the same radius, referred to as the micropore radius Rym; 

c) next molecules can be adsorbed in a direct neighbourhood of the site 
occupied by a molecule, provided that there is a space large enough, and 
the consecutive molecules are in direct contact with previously adsorbed 
ones (it may be treated as a multilayer adsorption of BET type). 

The assumption (a) is valid if a typical micropore wall distance is of range of 

double diameter of small molecule sorbates (ca 12:10"'° m) that is fulfilled for 

hard coals (refs [1-3]). It allows us to assume that the volume of the first layer 
adsorption space V,m is a property of the microporous structure, which does not 
depend on the sorbate molecule size (like Vasu for submicropores). Therefore, 
by virtue of the assumption (b), the surface of submicropores may be roughly 
evaluated, if Vm and R,, are given. In case of hard coals the radius of 
micropores R,,, is of less importance, as the main contribution to the coal 
surface is due to submicropores (refs [1-3]). However, for microporous sorbents 
this parameter may be of consideration. On the basis of our computations made 


Mathematical modelling of adsorptive-absorptive processes 45 


for an active carbon the value R,,,=8.5-10"° m (ca five times larger then water 
molecule radius) was found as appropriate [19]. 

The equation (8) describes also the pure adsorption (a=A; @,=0; w,=0), but 
under assumption that only Langmuir type monolayer adsorption occurs. To 
handle possible multilayer adsorption in microporous sorbents — see assumption 
(c ) — we assumed ([16, 18]) that there is a spectrum of adsorption sites in 
micropores, each being capable to contain the number of 1, 2, .. k, ... ce sorbate 
layers. Hence, creation of stack-like clusters (as in BET theory) is possible, but 
only for the limited number of k layers. In the paper [20] this concept was 
developed by admitting more compact clusters, in which the n-th layer of the 
clusters started from k-th type sites (containing m,,, molecules) provides mpin: Bin 
sites for the next layer (ın 21). In effect, we derived the following general 
equation describing the local adsorption equilibrium for the k-th type sites 
(constrained nonBET adsorption): 


* B? 
TB, * M okn-1 + (1+B;, )-m,., sate “Mons, = 95(14) 
kn+1 
where: n=1,...k; Mpo = myx (Mp Stands for the number of k-th type sites); 


Moyxk+1 =0, and 


Pot | 
A l E e ] 1-0-1], 
B= 8,,:|1-——-:- = ;  B,, =exp(Q,,/RT) (15) 
wr Mon, 

In the model (14, 15) proper representation of changes in configurational 
entropy is stressed [20], although side energetic interactions may be also taken 
ino account in the quantity Q'in. The set of eq.(14) has to be solved numerically 
(eg. by direct iterations) for specified sites (x). 

This approach seems to be useful for examination of adsorbents of regular 
porous structures (like zeolites), for which parameters mj, By, and Q",, can be 
specified. Thus, it offers an alternative way (not exploited yet) to the concepts 
presented in refs [22-24] for analysis of such adsorption systems. 

What concerns sorbents of irregular porous structures (like glassy polymers 
or hard coals), which are of our main interest, the model (14, 15) is not 
effective because of large number of parameters involved in, faced with highly 
uncertain parametrization of the sobent. Hence, in refs.[16, 18, 20] we assumed 
that the number of sites for consecutive k- layers adsorption may be roughly 
expressed by the following exponential function: 


mu= Mya-(1-0)-087 (16) 


where œ is an empirical parameter, mą- the number of sites for the first layer 
adsorption. 
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The above assumptions and certain simplifications (see refs [18, 20] bring 

the model (14, 15) to the following concise formula, referred to as LgBET 
equation (see ref.[20]): 
my (l-a) D-a- hen -p-0-a)] 
m, B+ B +11" -0-8 —a- BT") 
where: B=f,, = const for n=1, ..k; k=1,.. œ; B*, = B,/Bc; B,=B,,=const, 
B=B,,=const for n>1 (see eq.(15) with Q",, =Q, and Q",,= Qc =const for n>1, 
respectively), and: 
M'=1/B,; T” =m" -(-n")* (18) 
Nevertheless, we stated that the eq.(14) with B=1 is good enough to describe 
the constrained nonBET adsorption over relative pressures ranging to ca.0.85, 
and it is well applicable to typical adsorption systems with wide range of 
adsorption energies. In this pressure range effects of œ and Bc on m,a are hardly 
separable (isotherm is affected rather by o/B;). Hence, to complete MSM for 
sorption systems, we recommend application of the formula derived before in 
[16, 18], (referred to as LBET model), assuming that B=1 and BÆ1 (notice that 
B-=1 is commonly accepted in BET theory — see ref.[25]). The model LBET 
has the form: 


(17) 


m=] SOT g(a) 1 | TA 
i B, +1 1-a TI |B, (Q-mM+I 
B, =exp(Q} / RT) 20) 


where Q, '/RT is the sum of the all energetic terms in equation (8) applied to the 
A-th subsystem (ORZ). 


In spite of serious simplifications involved in the above formula, it seems to 
be advantageous, when compared to other constrained BET formulas (see 
e.g.[26]), as it represents clearly and relatively well effects of geometrical and 
energetic constraints on adsorption capacity. In turn, the model LgBET (egs.1 
7, 19) may be useful to explain singular properties of some adsorption systems. 
A convincing example is shown in Figure 2 [20], where LgBET and LBET 
formulas were applied to adsorption data measured and discussed by Choma & 
Jaroniec [27]. 

The same approach may be used to describe local adsorption isotherms on 
heterogeneous surfaces. In order to derive a heterogeneous adsorption formula 
based on LBET equation one may consider the parameter Z, in eq.(6) as 
uniformly distributed over the first layer adsorption sites. In the ref.[16] such 
a model was derived assuming that the distribution is the same for all sites 
(hLBET formula). It can be also done with a class of different distributions. 
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Let us take that the distribution function ranges from Z, (the same for all sites) 
to Zax, which is correlated with the number k of layers possible to be kept on the 


site: 
C 
Za =Z; Za =z, ($+a-0) (21) 


Ce (0,1) is a constant factor specifying the width of the energy spectrum. 


0.2 0.4 0.6 7 ars: 


Fig. 2. Sorption of N, on synthetic carbon Ambersorb 563 (circles - Choma & Jaroniec 
[27]) and isotherms obtained with LgBET (eq.17) and LBET equations [19], compared 
to BET and Langmuir (L) isotherms fitted to the initial adsorption data 


In effect, a specific distribution of the surface energy is obtained, depending 
on the site capacity distribution m,(k). Such a representation of the surface 
heterogeneity makes possible to integrate analytically the local LBET 
adsorption isotherms over the set of sites, provided that an approximation to the 
eq.(19) for local isotherms is taken. Possible approximating formulas are 
discussed in ref.[21]. Finally, the adsorption isotherm m,,/m),4 is expressed in 
the following form: 


Mpa. Mra) vo (Ee h | ee fat | eo 
m,, B, +1 i-a-T In(B,/B,) | By -(—1) +1 


where: By=B,(B,/B,/Bc)~"; B,=exp(U,/RT); Qa(Zs,U,) is the energy of 
adsorption on the most adsorptive first layer sites (k=1), U, — cohesion energy 
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of adsorbate, Qc — adsorption energy on the second and further layers, and 
W(a.,I1) is the following correcting function: 


a- (+0 -a-:1) (1-11) 


W (œ, 11) = 
eit a-(1+T-a@-M1)-(-M)-(1-a@)-d-a@-M)-[Ind-a@) -M- In - 


(23) 

The above model may be used to determine the number m,, of the first layer 

sites in porous structures, including the sites of very low adsortpivity. Such a 

parameter is a measure of a porous structure geometry, slightly affected by the 

adsorbate properties. Moreover, the model may be useful to examine energetic 
properties of adsorbents. 


3. MSM BASED INTERPRETATION OF SORPTION PROPERTIES 
OF MATERIALS 


The model (8) and (19 or 22) makes possible to calculate the theoretical 
isotherms of adsorption in submicropores, absorption in elastic regions of the 
sorbent matter, and the sorbent expansion ratio due to absorption and sorption 
in submicropores. It needs the system parameters to be assumed, i.e.: energetic 
parameters U,, ô; polarity factor C, and correcting function parameters Z4, Zz, 
Rg, (see eq.3); the number M,a of holes of particular size Rpa, parameters 
describing the structure of the sorbent matter: volume fraction and average size 
of crosslinked chains, quasi-crystalline domains, and non-crosslinked chains u,), 
Xer» Ud Xa Uez, Xe2} Xa» Xel » Xe2, COMpactness parameter Z,,y for quasi-crystalline 
domains (see eq.9), number of typical crosslinkages per crosslinked chain L,). 
Moreover, the molar volume V, of sorbate in the reference state has to be 
known. 

The pore structure may be specified in the form of pore size distribution 
Mpa=f(Ria) that enables us to compute the values for Mna. It was stated [14, 16], 
that the left truncated normal distribution may be applied to express distribution 
of submicropores, and one point distribution expresses adequately the amount 
of micropores. Therefore, data characterizing the pore structure are: volume of 
submicropores Vasu» their average radius R,q, minimal radius Ramin and 
dispersion Oga, volume of the first layer adsorption space in micropores Vym, and 
micropore capacity distribution parameter a (see eq. 16). 

Total volume V,,, of the “large” submicropores (determined by the normal 
distribution for R,2 R,) is viewed as providing a room only for the first layer 
adsorption, like the volume V,m of micropores — see assumption (a). Therefore, 
the number of sites m,, for the first layer adsorption is calculated as: 


Mra = (Vrs + Vim) Vs (24) 
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Let V, denote the total pore volume available for the monolayer adsorption: 
Vi = Visus + Vim (25) 


The volumetric adsorptivity of the micropores Vmi, viewed by the above models 
is 


V mic™ (Vnus + Vim) (1-a); (26) 
and the total volume of pores Vpores available for sorption, may be evaluated as: 
Vpores= Vasub- Vrs T V mic = Vhsub + Vus o ( 1-)+ V iml (1-0); (27) 


Let us notice that all the parameters listed above are of direct physical or 
geometrical meaning, hence their values provide information on more general 
properties of the sorbent and sorbate (molecular and porous structure, 
intermolecular interactions, cohesion forces etc.) that may be used for different 
purposes. 

The equations (8, 19) have to be solved together for all subsystems, with an 
assigned pressure P. In this way one obtains values of m,a(P/P,) for a sequence 
of P, i.e. the set of theoretical sorption isotherms for the specified subprocesses, 
in particular, the absorption and volumetric expansion isotherms. Then by 
summing up the sorption capacities m,, for all the subsystems, one may obtain 
the theoretical sorption isotherm m,(P/P,). 

The multiple sorption equation (8) has no analytical solution with respect to 
Ma. Hence, it should be solved numerically for a finite number of sorption 
subsystems specified according to the radius R, of pores. For this aim the range 
Ri€ [Rimin — Rp] may be divided into a number of sections with an average value 
Ra (in our calculations we specify 10 sections, i.e. we consider 11 subsystems 
including the pure absorption — R,o=0). The number of sorption sites Mpa in the 
a-th subsystem is calculated by dividing the volume of a-th subsystem holes 
(evaluated according to the normal distribution of R, within a-th section), by 
the volume V, of sorbate molecule. The last section, corresponding to the 
average radius equal to that of sorbate molecule (R,,=R,), is attributed to the 
pure adsorption subsystem (pointed with the subscript A). To this subsystem 
one also includes the first layer adsorption at remaining “large” submicropores 
and micropores. 

The model parameters are to be determined by fitting the theoretical 
isotherm to empirical sorption data for the examined system. In spite of a large 
number of the parameters the fitting is usually rather difficult. The quantities 
Uel, Xel » Ud Xa, Uez, Xe2; Xa aNd Zya describing the sorbent molecular structure and 
the solubility parameters 5, may be tightly bounded on the basis of literature 
data (see ref.[1]). MSM isotherms are strongly affected by the parameters of 
quasi-crystalline domains (uz Xa. Za), average size of submicropores (and 
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volume of micropores, and by energetic parameters (U,, 5.; Cp, Za Zs, Rp). 
Therefore, these parameters are well determinable. In particular sorption 
capacity is considerably affected (in energetic way) by the specific surface area 
of submicropores, which makes possible to detect a fraction of porous structure 
unavailable for sorbate molecules. Accurate description of elastic chains (u,;, 
Xen Len Xe2) is of less importance 

In general, the examination of a system does not provide precise values for 
the model parameters, but rather shows an area in the parameters space, in 
which the sorption data may be fitted. Nevertheless, this area may be highly 
tightened by examination of the sorbent with a number of different probing 
sorbates. Notice that the sorbent structure parameters (Visit, Vims Rhavs Orm Uel, Xer 
, Us Xd Uez, Xe2; Xd, Zwa ) and &, should not depend on the probing sorbate. The 
remaining ones may be bounded in such a way to force consistency of the 
examined systems parameters. Moreover, the results obtained with MSM may 
be confronted with data obtained in another way, especially, with sorption 
induced expansion of the sorbent (see refs.[12, 15]). 


empirical data 


pure absorption 


pure adsorption 


~~ first layer adsorption 


C, 
0) oinik op kaadaa eataa Oleh ak 


o 


Fig. 3. Sorption of CH3OH in a hard coal type 42 (92.4%C) (Wójcik [28]) and MSM 
isotherm calculated with eqs.(8, 19) 


Numerous examples of such multilateral studies of hard coals of different 
rank and type are shown in refs.[6, 7, 9, 10, 12, 15, 19, 28]. They give evidence 
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that sorption capacity of coal is affected by, both molecular and porous 
structure, and confirm significant role of submicropores in sorption processes. 
We showed that sorption induced swelling of coals may my well explained as 
an effect of absorption in elastic regions of coal and in submicropores. 

A representative example of MSM fitting results is given in Figures 3, 4 
(taken from [28]). 


empirical data f 
pure adsorption 


first layer adsorvtion--------~ 


Po 


Fig. 4. Sorption of HO in a hard coal type 42 (92.4%C) (Wójcik [28]) and MSM 
isotherm calculated with eqs.(8, 19) 


The figures show that sorption mechanisms of H,O and CH,;OH in coals are 
essentially different . The sorption of H,O is mainly by multilayer adsorption in 
microposes, although absorption is large enough to explain measurable coal 
swelling. In turn, CH;OH is sorbed mainly by penetration among elastic chains 
and filling of submicropores. 


4. CONCLUDING REMARKS 


The multiple sorption model may be viewed as a research instrument, 
supporting analysis of structures of elastic-porous materials. It yields 
a quantitative description of the material structure and intermolecular 
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interactions in sorption systems. Such a description may be used to answer the 
questions concerning the material structure, surface properties and mechanisms 
of sorption. The model needs more advanced computation tools, then 
commonly applied adsorption formulas (BET, DR). Nowadays it seems to be of 
no importance. 

Fitting of MSM to empirical data does not provide precise values for the 
model parameters, but shows an area, in which the sorption data may be 
explained with MSM mechanisms. Hence, the examination of a sorbent needs 
different probing sorbates to be used. Moreover, the results should be critically 
analyzed and confronted with data obtained in other ways. Nevertheless, the 
parameter determined in such a multilateral analysis may be used for prediction 
of the sorption system properties in pressure conditions unachievable in 
laboratories (e.g. analysis of hard coal sorption capacity and swelling in coal 
beds [15]). 

Apart from MSM we propose a set of formulas applicable to examination of 
adsorption systems. In particular LBET and LgBET models were found to be 
promising alternative for commonly accepted BET and DR equations. 
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